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ABSTRACT 


It is likely that supermassive black holes (M * 10^ to 10^ 

M ) exist in the nuclei of many quasars and galaxies. The collapse 

H 

which forms these holes and subsequent collisions between them 
should produce strong, broad-band bursts of gravitational waves: 
for a source of mass M at the Hubble distance of - 10^ liRht 
years, the dimensionless amplitude would be h ■ 2 * 10 ^ * 

(M/10 h M^) , and the duration of the burst would be r * (90 sec)* 

(M/IO^ M^) . Such bursts m*ght arrive at Earth as often as 50 
times per year — or as rarely as once each 300 years. The detection 
of such bursts may be possible within the next few years using 
dual-frequency doppler tracking of interplanetary spacecraft. 


Subjec t headings : galactic nuclei - gravitation - quasistellar 

sources or objects - relativity 
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I. SOURCES OF THE BURSTS 


When build Ing models for the cores of quasars, of strong radio sources, and 
of active galactic nuclei, theorists typically conclude that, whatever may be 
in the core, it is likely to generate one or more supermassive black holes 
(M * 10^ to lo ' 1 M ) in a time short compared to the age of the Universe. See, 
e.g., Lynden-Bell (1969), Wolfe and Rurbldge (1970), Spitzcr (1971), and Sasluw 
(1974, 1975). 

A typical scenario involves the collapse of supermassive stars, gas clouds, 
or star clusters to form holes, and perhaps subsequent collisions between holes. 
Each such collapse or collision involving masses M i 10^ M should oroduce a 
burst of gravitational radiation so strong that it might be observable at Earth. 
However, in the alternative scenario of a single hole which forms small (M << 

10^ M,.) and then grows bv gradual accretion, no strong burst of gravitational 
waves can be expected (Davis et aj.. 1971). 

There are so few quasars and galaxies in the Universe (~ 10^), and each 
one can produce so few strong gravitational-wave bursts (perhaps <; 10) that 
to detect several bursts per year, we must search for them over the entire volume 
of the observable Universe. Most such bursts will probably come from a redshift 
z = 2.5, when the Universe was about 2 billion years old, since quasar activity 
seems to have peaked very sharply at that time (Schmidt 1970, 1972). 

II. MEAN TIME BETWEEN BURSTS 

The mean time At between such gravitational-wave bursts at Earth can be 
expressed In terms of: (i) the Hubble expansion rate H () and deceleration parameter 
q f of the Universe today (which we assume globally Friedmann with zero cosmological 
constant); (il) the redshift z * 2.5 at which most of the hursts were generated; 
(lii) the number density n^ today of "centers" where bursts originated (each 
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center was probably a galactic nucleus); (lv) the mean number of bursts N generated 

In each center during Its active life; and (v) the speed of light c: 

1 [ 1 - q + q z - (1-q )(1 + 2q z) 1/2 j 

At - = , R ■ z ; (1) 

4»iR n Nc (H /c)q 'l ♦ z) 

o o ^o 

cf. equation (29 . 33)of Mlsner , Thorne, and Wheeler (1973) — cited henceforth as MTW. 
Because number densities of cosmological objects are usually computed using 

H ■ 100 km sec”' Mpc”' and q * 1 , we shall assume these parameteru throughout. 

o o 

The more modern values • 55 and q q « 0.5 (Sandage 1972; Gunn et^ al^. 1974) 

would give values of At smaller by a factor of 3. With H » 100, q ■ 1, 

° o o 

z * 2.5 the time between bursts is 

At • (0.06 years) (n Q /Mpc 3 ) * N * . (2) 

Various assumptions — each quite reasonable — lead to widely differing estimates 
of n Q , the density of burst centers. If nearly every quasar (radio-quiet or 
radio-loud) produced at least one black hole and accompanying burst, then n Q Is 
the density of dead quasars in the Universe: 


n * 

o 


h n llv„ (t) 

0 ‘life 



r 

\ Mpc 3 / 0 


10 5T di 


0.8 * 10 ' 4 \ 
Mpc 3 J 



Here »j jv> (t) is the comoving density of live quasars at time t with < -22 

as estimated by Schmidt (1970, 1972); tj lft is the mean lifetime of quasar activity 

in each center; and is the present age of the Universe. Schmidt's evolution 

law, n,, * lO^^o * ll*5t/t Q requires an upper limit t,., < t / 1 1 . 5 ; and 

live life o 

a lower limit might be t,^ e * (typical size, 100 kpc , of extended radio-emitting 
regions) /( typical velocity, U.t c, expected for expani ion of extended regions 
in models of radio sources) * 3 * 10^ yrs (cf. Longair, Ryle, and Scheuer 1973; 
Blnndford and Rees 1974). Thus 

1 x 10“ J Mpc’ 3 < n^ $ 0.2 Mpc 


if all quasars wete burst centers. 


(3) 


Suppose, alternatively, that all galaxies with violently active nuclei produce 

supermusslve holes and wave bursts during their violence. Presumably the violence 

seen In strong radio galaxies Is sufficient. This leads to an n^ roughly the 

same as for quasars (eq. |3]) (Schmidt 1966, 1972). 

Presumably the violence of a Seyfert galaxy also qualifies. The number density 

-4 - 3 

of Seyferts active today is * 2 * 10 Mpc (Huchra and Sargent 1973 adjusted 

from H^ ■ 75 n ■ 100). However, It is likely that Seyfert activity is transient, 

and that in the past a large fraction of all spiral galaxies with My < -19 experienced 

-3 

nuclear violence comparable to that in Seyferts (n^ - 0.01 Mpc ; Shapiro 1971). 

Arguments for such widespread violence In spirals include these: (i) Of all 

non-Virgo-cluster and non-local-group spiral galaxies with corrected recession 

velocities v < 1000 km/sec, * 10 per cent have small radio sources in their nuclei, 

and ' 5 per cent arc "optically active" in ways that could be remnants of earlier, 

more violent behavior (Ulrich 1974). (ii) The only big (My < -19) spirals close 

enough to be studied from the ground with optical resolution < 10 pc are our galaxy 

and M31 (Andromeda); and they both exhibit peculiar nuclear activity: (a) both 

show radial gas outflows suggestive of violent nuclear explosions in the past, 

8 •• 

perhaps 10 years ago (Munch 1960, Rubin and Ford 1969, Sanders and Prendergast 

1974, and references therein); (b) M31, as viewed by Stratoscope II (Light, 

Danielsen, and Schwarzschild 1974), has a nucleus so compact (M - 5 * 10^ M (>) ; scale 

height 0.5 pc) that all of its stars more massive than * 2 M 0 may by now have sunk 

to the center, producing violence (cf. Spitzer and Schull 1975); (c) the central 

3 pc of our own galactic nucleus contains a radio source with diameter $ 200 AU 

33 

and luminosity 3 10 erg/sec at v ' 10 GHz (Lo ej^ a_l. 1975), which coincides 

with a discrete near- infrared source (Becklin and Neugebauer 1975 source 16) that 

does not show CO absorption (Becklin, private communication) and that therefore 

could be a cluster of red dwarfs emitting . 3 10 J L and congregating around 

tota l “ 
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a supermassive black hole (Sanders and Lowlnger 1972); (d) one cannot rule out 

g 

our nucleus containing a black hole of mass a h large as 1 * 10 M Q (Sanders 
and Wr Ixon 1973) . 

If, as suggested by these observations, most spiral galaxies In the Universe 

have experienced enough nuclear violence In the past to produce supermassive holes 

.3 

and strong gravitational-wave bursts, n^ may be as large as 0.01 Mpc . If dwarf 
spirals > -19) also contain holes, but show little nuclear activity now 

-3 

because of a paucity of gas to accrete into the holes, then n Q could be 0.3 Mpc . 

On the other hand, if the only burst centers were the currently active Seyferts, 

-4 

then n * 2 * 10 . Tlius, 

0 

• A • 3 *3 

2 * 10 ' Mpc s n^ i 0.3 Mpc If galaxy nuclei are burst centers. (4) 

Note the striking resemblance of equation (4) to equation (3) — a resemblance 
which suggests to Lynden-Bell (1969) that most galaxies were quasars at one time 
and now contain supermassive holes. 

The average number of strong grav it at Iona 1 -wave bursts from each burst 
center is not likely to be large, perhaps 

1 < N s 10, (5) 

because the gravitational and collisional effects of each massive object are 
likely to inhibit the formation and evolution-to-collapse of other massive objects. 
However, our ignorance here, like our ignorance of n^, is enormous. 

By combining equations (2) - (4) we arrive at a range of "reasonable values" 
for the time between gravitational -wave bursts at Earth 

1 week i At i 300 years . ( 6 ) 

However, we must admit that hardly any strong bursts at all (At >> 300 years) is 
also reasonable; holes might usually form small and grow by gradual accretion. 


4 


III. EXPECTED CHARACTERISTICS OF THE HURSTS 

Consider the collapse of a massive star or star cluster to form a black 

hole, or the collision between tw< black holes with impact parameter of the 

order of their size. Let the total mass involved be M and let the total energy 

2 

radiated as gravitational waves be tMc . (A reasonable efficiency is e * 0.1) 
Most of the gravitational waves will come off in a burst with timescale 3/3 GM/r 
at the source (see, e.g.. Press 1971), and redshifted timescale at Earth 


r * (3/5 GM/c 3 ) (1+z) - (90 sec) (M/10 6 M y ) if z - 2.5. 


(7) 


The polarization-averaged dimensionless amplitude (h) 
waves at Earth can be computed fr> 


TT . TT . 1/2 - 

j h t j / 4) of 


the 


(e Me") - (total energy radiated) » ( ° | (4 hR 2 )t (l^z) 

\8hGt / 

( energy flux^ ^ ( surface area^ x (duration\ ^ /blue-\ 
at Earth J \around source/ lof burst) (shift/ ’ 


( 8 ) 


see equation (35.23) of MTV. Here R is as in equation (1). For * 100, 

q = 1, z ■ 2.5 this gives 
o 

(h) * (6/3 c) 1/2 CM/Rc 2 - 2 x 10 -17 (c/0. 1) 172 (M/10 6 M 0 ) . (9) 


Table 1 shows the expected (h) and t for sources uf different mass M. 

By observational studies of wave forms h(t) one may deduce much about the 

details of sources. For example: (i) If h(t) returns to its initial value after 

the burst, h,. . = h. ... ,, then the source involved only one object (e.g., 
final initial 7 J 

a single star or star cluster which collapsed to ferm a hole); but if h^^j t 
h . . . ,, then the source involved two or more discrete objects (e.g., two black 
holes that collided and coalesced). (ii) If part or all of the burst Is a rather 
monochromatic, sinusoidal signal, then the source may have been a rapidly 
rotating, relativistic star which went unstable by the "Chandrasekhar-Dedekind 
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bar-mode process" (Chandrasekhai 1970). (Ill) If part or all of the burst Involves 
a rather sinusoidal signal with a variety of harmonics present, then the source 
may have been a rapidly rotating, highly relativistic star which went unstable 
by the "Friedman (1975) ergotorold process". ( lv) If the burst Is entirely 
broad-band, except for a rapidly damped ringing at the end, the source may have 
been a star that encountc-ed a "dynamical instability" (see, e.g., Friedman and 
Schutz 1975), and collapsed to form a black hole. 

IV . DETECTION OF i'HF. BURSTS BY DOPPLER TRACK INC OF SPACECRAFT 

In the next few vears the best detector for the gravitational -wave bursts 
of Table 1 will probably be doppler tracking of one or more distant Interplanetary 
spacecraft: A tracking station on Earth transmits monochromat ic electromagnetic 

waves with frequency locked onto a highly stable clock (master oscillator). The 
spacecraft receives the waves, amplifies them, and transmits them back to Earth 
(l.c., it "transponds" them). The tracking antenna at Earth receives the 
transponded waves, compares them with the transmitted waves by counting the number 
of cycles and fractions thereof during an integration time tj nt * a f d derives a 
frequency shift from the comparison. (For further detail see, e.’ ., Anderson 
1973). The measured frequency shift is due primarily to the ve city of the 
spacecraft relative to Earth (doppler shift), but it also contains contributions 
("noise") from frequency fluctuations and drift of the master oscillator, from 
time-varying dispersion of the electromagnetic waves in the Earth's ionosphere 
and troposphere and in the interplanetary medium, and from passing gravitational 
waves . 

The gravitational -wave contribution as computed by Estibrook and Wahlquist 
(1975) can be described as follows. Orient the spatial axes of a Euclidean 
coordinate system so the gravitational waves travel in the z direction, and the 
Earth-spacecraft line of sight has length £. and lies in the x-z plane at an angle 
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ft relative to the z axis. Then the component of the gravitational waves which 

influences the doppler signal is h TT ; h(t - z). (For the meaning of h 1 ‘ see 

XX XX 

chapter 35 of MTV/. ) The influence of h on the doppler signal is given by 


Av 

v 


( 1 “ C08 q \ l _ . (• 1 + cc 

- { — i — ) h « - 6 S * — r 


(10) 


(Estabrook and Wahlqulst 1975). Here Av/v is the frequency shift of a specific 

short piece of the electromagnetic wave train — a piece received at Earth at 

time t ; and h , h._, b D are the values of h encountered by that same piece of 
K fc» l K 

wave train at its moments cf emission, transponding, and reception: 


hp. = h(t R - 20 , h^ = h (t R - l[l + cos ej), h R * h(t R ). (11) 

TVo features cf equation (10) are important: Fijrsj, the magnitude of the 

grav it at iona 1-wave- induced frequency shift is Av/v * h; second, (as Estabrook 
and Wahlqulst emphasize) the gr iv i tat Iona 1 wave-form is repeated three times in 
the Doppler signal, with relative amplitudes and time separations that depend 
on only one unknown parameter: the angle 0. This second feature may allow one 
to extract gravitat ional-wave effects in the presence of much larger doppler 
noise, and may also allow one to determine the angle 6. By simultaneous tracking 
of several spacecraft one could dig even deeper into the noise, and/or one could 
determine the 2-d i mens it lal direction of the source. 

When using single-frequency radio waves ("S-band", X * 14 cm), the NASA 

deep-space network can track spacecraft under favorable conditions (no planet 
or sun or high-density plasma cloud in the radio beam) with rms doppler residuals 
of 1 

*0ur remarks here and below about NASA doppler capabilities and possible improve- 
ments are b; ,<d largely on discussions with JPL personnel, including John D. 
Anderson, Richard Davies, Frank B. Estabrook, Richard Sydnor, Oldwig H. Vonroos, 
and Briant Winn. 
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i 7 min. 


( 12 ) 


<5(Av/v) * (3 * lO -13 ) for t. 

I n t 

where is the integration time. This precision is far short of that required 

for detection of the wave burst* shown in Table . (6(Av/v) » h * 10 to 10 
However — and this is the key point — there has not been much motivation for 
improvements in doppler tracking in recent years. Given the motivation, c ne 
could improve on the present precision by several orders of magnitude. 

The following improvements, largely motivated by non-doppler considerat ions, 
are underway or could be iitplc nented with small effort and cost: (i) conversion 
from s ngle-f requency tracking to dual-frequency ("S-band" from Earth to space- 
craft ("uplink"]; "S-band" and "X-band" (A * 3.6 cm| from spacecraft to Earth 

("downlink"!; conversion recently completed); (ii) replacement of rubidium clocks 
by hydrogen- maser ciocks as the master oscillators (replacement to be complete 

by summer 1976); and (ill) replacement of the present "doppler extractor," which 
-2 -3 

is accurate to 10 " cycles, by a new extractor accurate to 10 cycles (replacement 

not now planned, but easily doable and defin'tcly needed for the doppler system 

to take full advantage of the new master oscillators). With these improvements 

the system residuals could drop as low as 6(Av/v) * 2 * 10 ^ for t. > 15 nin — 

lnt 

though dispersion of the signal, even with dual -f requency monitoring, may make 

the overall residuals somewhat worse than this. By using the triplet structure 

of the gravitational-wave signal to dig Into the noise, and by transponding off 

spacecraft more distant than about 10 A.U., one might be able to detect 

-16 

waves with h * 1 * 10 and t - 70 min (case 11 of Table 1). NASA could 
try this in late 1976 with Pioneer 10 and 11 (the two spacecraft which flew by 
Jupiter in 1973-74); though such an experiment will be degraded by the lack of 
dual-frequency capabilities in the Pioneer spacecraft. NASA could also try with 
the Mariner Jupiter-Saturn mission (1977 launch; Jupiter flyby in 1979). 
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Further doppler improvement* might produce a capability for detecting 
case-ill waves (h * 1 "10 **’, t * 7 rain) and conceivably case-IV by late in 

this decade or early in the next. (For such waves one can use any spacecraft 
more distant than 1 A.U. — i.e., almost any Interplanetary spacecraft at all). 

Among the possible Improvements are: (1) X-band tracking on uplink as well as 
downlink (work on this begins in 1978); (ii) replacement of hydrogen -maser master 
oscillators by superconducting-cavity stabilized oscillators, which in 1974 had 
a stability A(Av/v) • 6 * 10 for > 10 sec (Stein and Turneaure 1975), 

which will probably reach 6(A v/v) * 1 * 10 lf> in 1976 (Turneaure, private communi- 
cation), and which — using niobium-coated sapphire cavities — may ultimately achieve 
stabilities much higher than l * 10 ' (Braginsky and Manukin 1974); (iii) the 
design and construction of new tracking-station electronics (doppler extractors, 
receivers, etc.) to go along with the improved master oscillators. 

Buffeting of the spacecraft by f’uctuating solar radiation pressure and 
solar wind will hide the gravitat iona 1 -wave signals at some level (perhaps below 
case III?) despite their unique triplet structure. If necessary, one can reduce 
buffeting by using a drag-free ("conscience-guided") spacecraft (Kundt 1974) or 
■i twin-probe spacecraft (Bertotti and Columho 1972). The TRIAD drag-free space- 
craft, which has already been developed, constructed, and flown once (Johns Hopkins 

-9 2 

and Stanford 1974) is sufficiently drag free (a < 5 * 10 cm/sec ) to reduce 
time-averaged buffeting noise to 6(Av/v) * 1 * 10 — a n °l 8e level 

comparable to the gravitational-wave signals of Table 1. However, the TRIAD 
"dead band" (maximum ^ at displacement between firings of correction Jets) was 
l millimeter — too large for gravitational-wave detection. A dead band of 0.1 p 
would be needed to c ompletely prevent buffeting at the level of case-IV wave bursts. 
A twin-probe spacecraft is not subject to such dead-band limitations. 
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Time-varying dispersion cannot be completely compensated for by dual- 
frequency ranging. One might ultimately remove dispersion in the Earth's atmosphere 
and near the Earth, as well as all noise sources at the tracking station (e.g., 
fluctuations In antenna dimensions), by tracking two spacecraft at different 
distances (e.g., 1 A.U. and 10 A.U.) but in the same line of sight, and by differ- 
encing their doppler signals (Davies 1974). And someday one might remove dispersion 
altogether by tracking at optical frequencies. 

Ore might he ,ole to search for case-11 and case-III bursts without the 
enormous expense of special-purpose spacecraft (drag-free, twin-probe, or two- 
in-1 ine-of-slght) . Detailed studies are needed to determine this. 

Simultaneous with any senrch for such grav Itat ional-wave bursts one should 
also watch for associated electromagnetic events (radio, infrared, optical, or 
X-ray) and for possible redshift changes due to mass outflow in gravitational 
waves. ’Ilie electromagnetic flux at Earth and apparent bolometric magnitude would be 


r - (o.3 -ip-) 


; E /E 
cm g w 

i i /t 
\ cm gw t 


m 


BOL 


-10.2 - 2.5 loSjol^T^ 

em gw 


(13) 


The duration of the electromagnetic event, t , should be much longer than that 

Ml 

of the gravitational-wave burst, t , because of absorption and reemission by 

gw 

gas surrounding the hole. The electromagnetic energy, E , might be far less 

CTT1 

than the gravitational energy, E — or they might be comparable. 

gw 

Even if the gravitational waves of Table 1 do not exist, doppler tracking 

-19 

should ultimately reveal waves at h * 10 due to (i) close binary systems in 

our galaxy; (ii) the collapse of normal stars (M < 100 M_) to form black holes in 

nearby galaxies (distance < 10 Mpc); (iii) the capture of globular clusters, 

3 4 

containing black holes of 10 - 10 M , by galactic nuclei at a distance -1000 

Mpc, followed by collisions and coalescences of the holes in the nuclei (cf. 
Bahcall and Ostriker 1975; Silk and Arons 1975; Tremaine, Ost r iker, and Spitzer 
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1975); And/or near encounters between compact stars In a possible dense star 

cluster In the nucleus of our own galaxy (7.el'dovlch and Polnarev 1974). In the 

case of close binary systems the periodicity of the signal can be used to enhance 

-18 

the uccuracy of the experiment to an effective level of h * 10 (Davies 1974, 
Estabrook and Uahlquist 1975). 

V. SUGGESTIONS FOR FURTHER STUDY 

In view of the above discussion, we suggest the following: (1) Tliat space- 
research organizations (NASA, USSR Academy of Sciences, ESRO) initiate detailed 
studies of the expected limits which various doppler-tracking configurations can 
place on gravitational waves of the type shown in Table 1. (2) That, if those 

studies are as optimisMc as our cursory analysis, high priority be given to 
doppler-trai king se mm for gravitational -wave bursts. (3) That astronomers 
and astri pnyslclsts try to firm up our estimates of the density n^ of centers 
where strong wave bursts were generated, and of the number N of bursts per center. 
(4) That relativity theorists develop a catalog of gravitational wave forms to 
be expected from the various objects and events that might occur in the nuclei 
of galaxies and quasars. (5) 'Hint experimenters try to develop alternative 
detection schemes for the waves shown in Table 1. 

For helpful discussions or correspondence we thank the JPL personnel listed 
in footnote 1; also Eric Becklin, Peter Bender, Bruno Bertotti, Daniel de Bra, 

Peter Goldreich, Wolfgang Kundt, Charles Misner, Gerry Neugebauer, Jerry Ostriker, 
Wallace Sargent, Bill Saslaw, Maarten Schmidt, Lyman Spitzer, John Turneaure, 
and Tot./ Tyson. One of us (KST) thanks the Aspen Center for Physics for hospitality 
during part of this research. 
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TABLE 1 


CHARACTERISTICS OF CRAVI TATIONAL-WAVE BURSTS FROM BLACK-HOLE EVENTS 
INVOLVING A MASS M AT A REDSHIFT z - 2.5 


Case I Case II Case III Case IV 


Mass , 

M 

5 * 10 8 H 

H 

5 * 10 7 M„ 

f) 

5 * 10 6 M ft 

5 x 10 

Wave « 

•mplitudc, h 

, . 10-“ 

. . 

1 - io - 16 

1 x 10 

Burst 

duration, t 

700 min. 

70 n.in. 

7 min. 

40 rjc 
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